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bJQ In this paper, we want to study, in a self consistent way, some nonlinear and collective transverse effects due to the 
^3 interaction of a relativistic electron/positron beam with a magnetoactive plasma in overdense regime, i.e., uq » rib 
(no and rib being the unperturbed plasma and beam densities, respectively). To this end, in the long beam limit, 
we consider the self-interaction that is produced when a relativistic charged particle beam is travelling in the plasma 
^ exciting large amplitude plasma waves, namely plasma wake field (PWF) excitatioiP'. 
I— I Hereafter, we refer to our papei^ appearing in this proceedings, as well. According to this paper, we assume that 

the plasma is collisionless and cold, with ions at rest forming a uniform background of positive charge. Furthermore, 
Qh a strong constant and uniform external magnetic field is assumed to be acting along the z-axis, Bq = BqGz. We also 
^ assume that the electron/positron beam is initially travelling along the direction of the magnetic field with a velocity 
C V5 = (3cez — 1). We consider a fiuid model, consisting of Lorentz- Maxwell system of equations for the beam-plasma 
^ system. From the perturbed Lorentz-Maxwell system, we obtain an equation that governs the evolution of the plasma 
wake potential driven by the charged particle beam density. On the other hand, by ignoring the longitudinal beam 
dynamics, we write the equation that governs the spatio temporal evolution of the charged particle beam given by 
the thermal wave model (TWM)^ ^ (For details, see Ref.^). These pair of equations can be cast as a quantumlike 
Q Zakharov system of equations which governs the self consistent spatio temporal evolution of the PWF self-interaction 
of the electron/positron beam, viz., 

^ *e|| = -^Jvivl/-^i.(rxxVx)* + f/»(rx,6* + ^-?^^i*, (1) 
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\^ where uouh is the electron upper hybrid frequency and ^ = ^(r^,^) is the beam wave function (BWF), so that its squared 
modulus is proportional to the beam density, i.e., y9b(r^,<f) = {N/cjz) |^(r^,<J)|^, where N and Gz are the total number of 
particles and the beam length, respectively, K = {(jJc/'2joc)'^ = (^Bo/2mo7oc^) = (A:c/2)^, V5_ is the transverse part of the 
' I gradient operator, Uw = {Aiz — /mo7oc^ is the dimensionless wake potential with Aiz = Aiz{r±,^) and = (/)i(r^,^) 
QQ the longitudinal vector potential perturbation and electric potential perturbation, respectively. Here, is the transverse 
position vector, ^ = z — /Set ^ z — ct plays the role of time-like variable, mo and 70 are the electron/positron rest mass 
T-H and the unperturbed relativistic gamma factor of the single particle of the electron/positron beam, respectively. In cylindrical 
^— H coordinates, r±, cp, ^, we look for a solution of the Zakharov-like system of the form ^(r^, cp, ^) = exp [im (cp — kc^/2)] il^m{r±, ^) 
^ with m integer, taking the limiting case |v5_| ^ ujpe/c^^uH- L^t us define the transverse beam size in the form of r.m.s., i.e., 
crm(0 — r\\iljm\^r±dr±. Under the above assumptions and definitions, from the Zakharov-like system, we easily obtain 

the following 2D Gross- Pitaevskii-type equation, viz., 

^ .^V^m 1 d ( d^m\ r I / ,2 / ^ 2 ^ 

^^ = -2^a^r^a^j-^-l^-l ^-+^2^^^^+2^J^-' 

where we have introduced the following dimensionless quantities: ^ ^ ^//3o, r± r±/ao, ipm \/ TvmlaQ ipm, Kb = Kcfq/c^ ^ 
5m — nbo-o/no'yoe^ml, ao and e being the initial transverse beam spot size and the transverse emittance, respectively. 
We use the virial equation associated with eq. ([s]) to get the envelope equation d^{cF^^Y /d^^ + 4i^b(cr^)^ = 4:Am, where 
cr^ VmJ.alP^/ao and Am = |(m + 1)!(1 + Kt) - (5m(2m)!2~^*^'^+^^ is constant of motion. Note that, from the envelope 
equation, the matching condition for the equilibrium transverse beam spot size, cr^, is Kh{(J^qf' — Am- 

A preleminary numerical analysis has been carried out by solving eq. (Is]) assuming the initial normalized BWF (density profile) 
as ipm{r±, 0) = exp (r±/2^ . The spatio-temporal evolution of \ipm\ has been investigated for different values of m, Kb, and 
Sm, at ^ = 0, 0.25T, 0.5T, 0.75T, T, where T = tt/^/K^. For both m = and m = 1, when the matching condition of the 
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^ = = 0.9069 e, - 1.8138 ^ = 2.7207 e, = 3.6276 




^' = ^' = 0.785398 ^' = 1.5708 ^' = 2.35619 ^ = 3.14159 




FIG. 1: 3D plots of |'0m|^ as function of x/cr and yja for different values of the dimensionless time ^ for m = 0: K\y = 0.25, 
8m 0.5, 0.5 (first row); i^b 0.75, 8m 0.5, 0.75 (second row); K^, 1.0, 1.5, Am 0.625 (third row). 



^"= 1.5708 ^=3.14159 4^= 4.71239 ^ = 6.28319 




4=0 4 = 0.811156 ^' = 1.62231 ^ = 2.43347 ^' = 3.24462 




< = ^ = 0.641275 ^ = 1.28255 ^" = 1.92382 ^ = 2.5651 




FIG. 2: 3D plots of \\\)m^ as function of x/cr and yja for different values of the dimensionless time ^ for m = 1: = 0.25, 
8m = 0.5, Am = 1.1875 (first row); i^b = 0.9375, 8m = 0.5, = 1.1875 (second row); i^b = 1.5, 8m = 3.5, Am = 2.0625 (third 
row). 



envelope equation is satisfied, the profile is practically unchanged (see the second row of Figures 1, 2 and 3, respectively). This 
predicts the existence of nonlinear coherent states (sometimes called 2D solitons). Furthermore, due to the strong nonlinearity, 
the effect of beam halo has been observed for m = 2, as displayed by both 3D and density plots in Figure 4. Due to the 
interplay between the strong transverse effects of the plasma wake field (collective and nonlinear effects) and the magnetic field, 
envelope oscillations with weak and strong focusing and defocusing have been observed for m = (see the first and the third 
row of Figure 1), m = 1 (see the first and the third row of Figures 2 and 3, respectively) and m — 2 (see first and second row 
of Figure 4) . Finally, the existence of vortices (effect of the orbital angular momentum due to the external magnetic field) are 
clearly shown in Figures 2,3 and 4, respectively. 

The present investigation seems to be useful for the plasma-based focusing schemes to be employed in the final focusing stages 
of linear colliders as well as for manipulating relativistic electron/positron beams that suggests the new fields of nonlinear and 
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FIG. 4: 3D plots (first row) and density plots (second row) of |'0m|^ as function of xja and yja for different values of the 
dimensionless time ^ for m = 2: K\) — 1.5, bm — 5.0, Am — 5.625. 



collective singular electron optics. 
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